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Abstract: Background oriented schlieren (BOS) visualization technique is 
examined by means of optical geometry. Two most important results are 
the calculation of the sensitivity and spatial resolution of a BOS system, 

y^ \ which allows for the determination of the experiment design space. A set of 

relations that characterize the performance of a BOS measurement is given, 

CIh' with emphasis on the design of background pattern and spatial placement 

q 
o 

>»: 

^ ■ 1 Introduction 

Background oriented schlieren (BOS) is a flow visualization technique that 

^ I carries the promise of quantitative measurements, as introduced by Meier 

^ ; (2002) and Richard and Raffel (2001). The physical principle on which BOS 

0^ I relies is similar to that of the formation of mirages: a scene viewed through 

• i a hot plume of air will appear distorted. The distortion is caused by the 

f^ , bending of the light rays from the straight line propagation due to inhomo- 

^' geneities in the medium, and it is described by the functional dependence of 

the refractive index on space, n(r) (Hecht 2002). One can consider BOS as 

an improvement and/or simplification of defocused grid schlieren (Vasil'ev 

5-H ■ 1971), white light speckle photography (Giglio et al 1980), moire deflectom- 

" " etry (Merzkirch 2007), and similar techniques. A large group of researchers 
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refer to BOS as synthetic schlieren, with the main difference between the two 
being that BOS is more general in terms of the background, while synthetic 
schlieren employs an optimized synthetically generated pattern as a back- 
ground (Dalziel et al 1998). There are two fundamental differences between 
BOS and standard (Toepler) schlieren techniques: in the latter, one records 
the variation of illumination in the image plane, and also one intervenes on 
the light path, usually with a knife-edge, modifying it in addition to the 
modifications by the flow. 

In BOS, a pattern is recorded at two different instances, once without 
and once with the flow that is to be investigated, giving the reference and the 
measurement images, respectively. The flow causes inhomogeneities in the 
medium, which will influence the refractive index. Thus, when the reference 
image is recorded, the refractive index of the medium is considered to be 
constant. On the other hand, when the the measurement image is recorded, 
the refractive index is a function of space in the domain defined as the test 
section. Interrelationship between the refractive index and flow properties 
- mainly density, as expressed by the Lorentz-Lorenz equation - allows for 
determination of the state of the fluid. ^ Based on this, BOS has found 
many applications, as described in papers by Venkatakrishnan and Meier 
(2004), Kinder et al (2007), Moisy et al (2009), Mizukaki (2010), Kirmse et 
al (2011), Sourgen et al (2012), Glazyrin et al (2012), and Vinnichenko et 
al (2012), to name a few. 

BOS technique consists of two independent stages: image recording and 
image evaluation. The first deals with the problems of designing the pat- 
tern that is to be imaged, placement of the components in the optical setup, 
and determining appropriate conditions for recording the images. The im- 
age evaluation part deals with the problems of comparing reference and 
measurement images, and data extractions from those comparisons. The 
dominant method for data extraction in BOS relies on cross-correlation al- 
gorithms that have been developed for speckle photography and particle 



■''Lorentz-Lorenz equation relates polarizability of an assembly of molecules to the re- 
fractive index. For gases, this relation can be approximated by the Gladstone-Dale equa- 



tion (Born and Wolf 2005). 
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Figure 1: Optical setup for background oriented schlieren. 



image velocimetry. Alternative methods include optical flow (Atcheson et 
al 2009), laser interferometric computed tomography (Ota et al 2011), and 
recently some developments on single pixel correlation are being proposed 
(Kahler et al 2012). 

There have been several publications, including those by Elsinga et al 
(2004), Goldhahn and Seume (2007), Yevtikhiyeva et al (2009), Ambrosini 
et al (2012), and Gojani and Obayashi (2012), which assessed the sensitivity 
and accuracy of BOS technique, mainly from the data processing and image 
evaluation point of view, because these introduce the dominant measurement 
uncertainties. Nevertheless, in order to determine the performance of a BOS 
system, it is necessary to investigate the measurement uncertainties arising 
from the optical setup, as well. This becomes more important in the cases 
when experimentation is carried with instruments of modest specifications, 
such as the case of high-speed imaging, limited by the relatively low pixel 
count of the image sensor. 



2 Sensitivity 

BOS measurement is able to quantify ray deviation, but only integrated 
through the line of sight. This is illustrated in figure 1, which shows the 
meridional plane of a typical BOS setup consisting of a structured back- 
ground, the test section under investigation, (also referred to as the phase 
object, transfer channel, density field, etc.), the objective lens with focal 
length / (focal points F), the image recording sensor, and the respective 
distances. When the reference image is being recorded, a feature from the 
background located at the point B is imaged in the point /. Introduction 
of the fluid flow with variable refractive index n(r) will deflect the beam for 
an angle e, thus the imaged point now will be shifted for Ajj to the point 
J. A necessary condition for this effect to be observed is that the pattern 
shift Aij be larger than the pixel linear dimension ipx- Thus, the detection 
limit is determined based on the relation^ 

A consequence of the change of the ray direction is that the ray will 
exit the test section at a different point, as compared to the ray during the 
recording of the reference image. But, this difference, namely 5e in figure 1, 
is negligible for as long as 6e, which is the distance between the entry points 
of these rays in the lens, is much larger. Since 6e/Si ^ w/{so — st), this 
condition is satisfied if 

" «1. (2) 



So - St 

This means that in the case of a BOS setup with dimensions much larger 
than the width of the test section, one can assume that the fluid flow only 
deflects the light ray, but does not displace it. Furthermore, similarly to the 
thin lens which is described by the principal plane, the entire test section 
can be approximated by a refractive plane. 



This equations holds for a ray of hght that images a point from the physical space 
into a point in image space. As it is known from image analysis, the determination of 
pixel shifts with subpixel resolution is possible, but only over an area. 



The sensitivity S" of a measurement setup is defined as the ratio of the 
detected change - for the case at hand, the pattern shift Ajj, - to the corre- 
sponding change in the parameter being measured, - angle of deflection e, - 
that is, 

S=^. (3) 

Let the distance of a background point B from the optical axis be r^ (in 
the object space), as shown in figure 1. The lens with diameter of aperture 
A forms the image of B at point / by collecting only rays emitted within 
the range of angles a, for which 

n - A/2 n + A/2 

< tana < , (4j 

where a is the angle of the direction of the ray with respect to the optical 
axis, and Sq is the distance between the principal plane (the plane of lens 
and aperture) and the background. In practice, A ~ 1 — 50 x 10~^ m and 
So ~ 0.1 — 100 m.' For the paraxial approximation to be applicable within 
an accuracy of 1% {a < 15° ^ 0.25 rad), the relationship between r^ and Sq 
is 

rb< 0.25 So, (5) 

which means that the dimensions of the field of view should not exceed the 
quarter of the dimensions of the BOS setup. 

In BOS, the thin lens relation 1// = 1/sj -|- 1/so holds, with the lens 
being focused on the background. The image / is at a distance rj = Mr^ 
from the optical axis (in the image space), where M = Si/so is the optical 
magnification of the system. Here, since BOS is basically a photographic 
technique of real and inverted images, the magnification is taken to be pos- 
itive, as is customary in photography. In the case of the defiected ray, the 
image point J is at a distance rj = M[rt — {a + £)st] = ri — M e st, where 
rt = ri, + a St is the distance of the deflecting point T from the optical axis, 
St is the distance between the background and the test section (see figure 



^There are examples of using BOS technique for visualizing open air blast waves with 
So ~ 400 m and St ~ 100 m (Mizukaki et al 2012). 



1). 

Since Ajj = \ri — rj\ = e M st, the sensitivity can be expressed as 

S = Mst = ^^. (6) 

•So 

3 Spatial resolution 

In BOS measurements, reference and measurement images are taken under 
the same Hghtning conditions of the field of view, which consists of a pattern 
with highly varying brightness from point to point. Nevertheless, irrespec- 
tive of the radiant exitance from a background point, it is the irradiance on 
the pixel area - resulting in the image of the said point - that is actually 
used for image evaluation and quantitative measurements. Thus, for a linear 
proportionality between exitance and irradiance, it is the relative difference 
in irradiance from point to point in the image space that determines the 
pattern. 

BOS imaging is done with charge-coupled device (CCD) or complemen- 
tary metal-oxide semiconductor (CMOS) image sensors, which digitize the 
output of the light falling onto the area of a pixel. Digitalization of the 
image of the field of view leads to discretization in two domains: (i) in the 
space domain, an operation referred to as sampling, and (ii) in the domain 
of intensity values, an operation referred to as quantization. Quantization 
is influenced by several noise sources in the image sensor, but the following 
discussion will assume that the variation of irradiance due to imaged pattern 
is much larger than the combined noise levels (dark frame and individual 
pixel response). In a further simplification, it is assumed that the image is of 
a binary type; hence, when the (relatively) high irradiation covers more than 
half of the pixel, the pixel records light (pixel grayscale value 1), otherwise 
it does not (pixel grayscale value 0). 

Similarly to the previous section, let the image of the background point 
B, i.e. point /, be at a distance rj from the optical axis, falling on the i — th 
pixel (i = 1, 2, . . . , A^, where N is the number of pixels of the sensor in one 
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Figure 2: Cone of light creating the image in BOS. 



direction) . Because of samphng, all the other points that fall in the segment 

(7) 



^ 7T ) ^P'X S: fi 'Si I ^ "•" Q / ^px 



are indistinguishable (for analysis in 2D, the factor multiplying the length 
of a pixel should be iibl/v2, but this correction is insignificant in our case). 
Furthermore, let A^ = kipx be the diameter of the circle of confusion (CoC), 
with k a small integer (let us say, A; = 1, 2, ... , A^/100 pixels). Based on the 
properties of CoC, only irradiance variations that are at a distance larger 
than this diameter can be resolved, therefore all points within the segment 
Tj lb At/2 are indistinguishable. Then, condition (7) is modified to 



fe + 1 



< 



< 



''PX 



k + 1 



(8) 



where [x~\ designates the nearest integer function of x. 

The cone of light collected by the aperture with diameter A, as it is 
illustrated by the shaded area in figure 2, does not contain independent 
information along its cross-section, because image sensor records the average 
irradiance contained in it, thus it appears uniform. Backprojection of the 



segment defined by equation (8) through the lens into the background plane 
determines the smallest distance that can be resolved, 6{y, which has the 
value 

(9) 



rC -\- 1 tpx 



M 2 ' 

and this defines the diameter of a dot in the background, i.e. the shortest 
distance between pattern variations that are resolved. In practice, the dot is 
made two or three times larger, so that image evaluation by cross-correlation 
can yield an optimized result. 

The spatial resolution of a BOS system at the deflection point T corre- 
sponds to the diameter of the projection of the cone of light on the refractive 
plane, and its value is 

M St At A st\ 

+ T7 1-— ' (10) 



M + 1 /# M \ s. 

where /# = f/A is the focal ratio (f-number) of the lens, and where st/so < 1 
holds. 



4 Discussion 

In BOS, the test section is placed between the background and the objective 
lens. This limits the values of st in the range < st < Sq, and consequently 
the sensitivity of the BOS setup is 

0<5<Si. (11) 

This results shows that the sensitivity of a BOS setup can be adjusted to any 
desired level, limited only by the focal length of the lens. In practice, though, 
sensitivity is adjusted by deciding on the spatial extent of the experimental 
setup (adjustment of st, limited by the laboratory space and the field of 
view) and the lenses and camera used (adjustment of M). Older publications 
that treated the sensitivity of BOS argued for the use of long focal length 
lenses for increased sensitivity. While this is correct, the present paper gives 
a more fundamental reason on why long focal length lenses are preferred: 



longer focal length lenses result in higher magnifications. Nevertheless, as 
long as high magnification of the field of view (test section) is achieved, even 
lenses with relatively short focal lengths can be used. 

Spatial resolution of a BOS setup, which is a characterization of the 
blurred out-of-focus imaging, is influenced by two terms composed of sev- 
eral parameters, as given in equation (10): it can be said that Aj and /^ 
are features of the camera and have a limited range of values, typically 5 - 
50 //m for At, and 5.6 - 32 for /^. For very low values of magnification, 
such as in the cases of BOS applied in outdoor measurements (natural back- 
grounds), the dominant term tends to be the second one. But, these cases 
also involve large distances between the background and the phase object, 
with St usually of the order of a few meters. Hence, both terms would be 
accountable. For relatively large magnifications (M > 0.1), the second term 
can be neglected, by virtue of having a very low value for A^. Furthermore, 
when considering error propagation, even a small uncertainty in sj, say 0.1% 
(1 mm in a meter), would yield a larger uncertainty in 5t than the value of 
Af/M itself. Thus, a good approximation of the spatial resolution in terms 
of experimentally measured parameters can be expressed as 

Si St ^^2) 



Si + So f# 

This result is to be interpreted as follows: the angle of deflection e at the 
point T is influenced by the variation of the refractive index n(r) along the 
entire cone of light with the vertex at B and base diameter 6t, i.e. r = rt±6t- 
Standard schlieren techniques use the phase object as the imaging object, for 
which St = 0. Hence, they have a superior spatial resolution as compared 
to BOS, because only the term At/M contributes to the uncertainty of 
measuring the angle of deflection. 

Equations (6) and (10) are the main equations that characterize a setup 
for a BOS measurement, and their graphical presentation is given in figure 
3. The surface plot for St is shown for unit values of At and /#, and where 
St « So- This figure presents the BOS experiment design space, meaning 
that once the acceptable range of values for resolution and sensitivity are 
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Figure 3: BOS experiment design space. 

decided, one can determine the distances between optical components when 
setting an experiment. 

This simple treatment of a BOS setup did not take into account several 
other factors that influence the measurement, such as lens aberrations, de- 
viations from paraxial approximation, sagittal rays, and alike. But, some 
of these influences can be cancelled in the image evaluation stage, under 
the condition that the recording of the reference and measurement images 
is done with a fixed position of the camera and the test section relative to 
each other. Also, the irradiation model was very simple and did not consider 
the varying response of the image sensor. In addition, considering the small 
size of a pixel in image sensor, the diffraction effects may play a crucial role 
in imaging, especially if lasers are used for illumination.^ Nevertheless, the 
obtained results are quite general and can be applied to several variations 
of BOS setups. An important result is that neither S or 6t depend on the 
angle a (admittedly, as long as condition (4) is satisfied). 



In order to eliminate diffraction effects, equation (9) with fc = may serve as a 
threshold for the coherence length of the light source. 
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5 Conclusions 

The present paper is a contribution in examining the performance of back- 
ground oriented schheren (BOS), and it supplements our previous work in 
assessing BOS (Gojani and Obayashi 2012). Several guiding formulas for 
setting a BOS system are presented: 

• detection limit is given by formula (1), 

• condition for no-displacement of ray at the exit of test section is given 
by relation (2), 

• dimensions of a field of view for a given spatial extent of a BOS setup 
is limited by relation (5), 

• system sensitivity can be estimated by equations (6), 

• the minimal diameter for a dot in the background is given by equation 
(9), and 

• a fairly good approximation of spatial resolution from measurement 
setup is calculated by equation (12). 

These formulas should be used as a starting point in setting a BOS 
system and determining the needed instrumentation, but they are not con- 
clusive for the latter, because in this analysis, the quantum efficiency of the 
image sensor and the temporal requirements of the experiment are not taken 
into account. The emphasis of this study, instead, was on the spatial setting 
of the instruments, so that desirable effects can be detected. 
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